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Abstract. This paper presents results obtained from Stokes I and V spectra of the B2Vp star 
sigma Ori E, observed by both the Narval and ESPaDOnS spectropolarimeters. Using Least- 
Squares Deconvolution, we investigate the longitudinal magnetic field at the current epoch, 
including period analysis exploiting current and historical data, a Ori E is the prototypical 
helium-strong star that has been shown to harbor a strong magnetic field, as well as a magneto- 
sphere, consisting of two clouds of plasma forced by magnetic and centrifugal forces to co-rotate 
with the star on its 1.19 day period. The Rigidly Rotating Magnetosphere (RRM) model of 
Townsend & Owocki (2005) approximately reproduces the observed variations in longitudinal 
field strength, photometric brightness, Ha emission, and various other observables. There are, 
however, small discrepancies between the observations and model in the photometric light curve, 
which we propose arise from inhomogeneous chemical abundances on the star's surface. Using 
Magnetic Doppler Imaging (MDI), future work will attempt to identify the contributions to the 
photometric variation due to abundance spots and due to circumstellar material. 
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1. Introduction 

The term "magnetosphere" was coined to describe "the region in the vicinity of the 
earth in which the Earth's magnetic field dominates all dynamical processes" (Gold 1959). 
The Earth's magnetic field is distorted by the solar wind, compressed on the side closest 
to the sun and pulled out on the opposite side into the magnetotail. As observations 
improved, it was shown that each planet with a magnetic field also possessed a similar 
magnetosphere to the Earth's, all created by the solar wind's effect on the magnetic field. 
Eventually, "magnetosphere" was expanded to include stellar environments where winds 
and magnetic fields interact, e.g. the solar win and magnetic field. 

Unlike planetary magnetospheres, stellar magnetospheres expand outward, as the wind 
is coming directly from the star itself. We can understand solar system magnetospheres 
in great detail, describing the complex plasma physics, as we are able to obtain direct 
observations by satellites. In the stellar case, we must use remote observations to infer 
physical properties (i.e., temperature, density), augmented by MagnetoHydroDynmaical 
(MHD) simulations. 

There is a distinctive difference between solar and massive-star magnetospheres. In 
lower mass stars like the sun, an a-Q dynamo combines a convective envelope with dif- 
ferential rotation to produce complex, small-scale, time variable field structures. Globally, 
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the field is quite weak. The solar wind is pressure-driven with a relatively weak mass loss 
rate (~ 10 -14 M /yr). More massive stars have mostly radiative envelopes. This major 
difference makes magnetic field generation by a dynamo difficult to explain, so the ques- 
tion of the origin of magnetic fields in massive stars remains open. Massive-star fields are 
also quite different from their solar counterparts. Their fields are generally steady, large 
scale, and sometimes very strong. Massive-star winds are radiatively driven with strong 
mass loss rates (on average ~ 10 -7 M /yr). Hereafter, this paper will focus solely on 
massive-star magnetospheres, specifically in the strong magnetic field case. 



2. Massive star magnetospheres 

As the wind interacts with the magnetic field, material from the stellar wind can 
be perturbed, channeled, torqued, and even confined. Massive star magnetospheres are 
generally structured and dynamic. Material is released in reconnect ion events, falls back 
onto the star and creates collisional shocks. Rotation complicates the physical picture, but 
it also allows for regular modulation of observables, through which the physical properties 
of these magnetospheres can be identified. The magnetospheric properties depend on the 
relative strengths of the wind and magnetic energy densities of the star. ud-Doula & 
Owocki (2002) define a wind confinement parameter to quantify this struggle between 
the wind strength and magnetic strength, 

D, it, 

V* = (2-1) 

For 77* < 1, the wind overwhelms the field, but for 77* > 1, the magnetic field is dominant 
near the star, channeling wind material into a magnetosphere. The Alfven radius, Ra, 
is where the magnetic and wind energy densities become equal. For a strong magnetic 
confinement (77* ^> 1) and a dipole field configuration, Ra ~ 77* 1 ^ 4 R* (ud-Doula & 
Owocki 2002). ud-Doula et al. (2008) also define a Kepler co-rotation radius (Rk = 
(vrot/vcrit)~ 2 ^ 3 R*)- If the Alfven radius is father out than the Kepler radius a region of 
confinement is created. Wind material is spun up by the magnetic field and supported 
by centrifugal forces. Above Ra, the field is too weak to confine material; below Rk, it 
is no longer supported by centrifugal forces and so falls back on the star. 



3. The RRM model and a Ori E 

In many magnetic Bp stars, the magnetic field is strong, rotation is fast, and 77* ^> 1. 
The field spins up and confines the channeled wind keeping it rigidly rotating well beyond 
the Kepler co-rotation radius, but also held down against the net outward centrifugal 
force relative to gravity. For this strong field limit, Townsend & Owocki (2005) developed 
an analytical model, the Rigidly Rotating Magnetosphere (RRM) model, wherein wind 
plasma settles at the local minima of the effective (gravitational + centrifugal) potential. 
The plasma then accumulates in the magnetosphere and is forced to co-rotate with the 
star. 

a Ori E (HD 37479) is a helium-strong B2Vp star with a 10 kG global dipole magnetic 
field. It is an oblique rotator with two plasma clouds rigidly rotating with the star (Land- 
street & Borra 1978) with rotational speeds greater than the surface speed of v sin i = 150 
km s _1 . Observations show modulation according to the 1.19 d rotation period in longi- 
tudinal magnetic field (Landstreet & Borra 1978), Ra emission (Walborn 1974), He line 
strength (Pedersen & Thomsen 1974), photometry (Hesser et al. 1976), UV line strength 
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Figure 1. Figure from Townsend et al. (2005). Observed (diamonds) and modeled (solid line) 
Stromgren u-band light curves of a Ori E, phased on the star's 1.19 day rotation period. 

(Smith & Groote 2001), 6 cm radio emission (Leone & Umana 1993), and linear polariza- 
tion (Kemp & Herman 1977). a Ori E also has anomalous surface abundances of He and 
Si (Reiners et al. 2000). The RRM model reproduces the photometric, magnetic field, 
and Ha observations relatively well. It only accounts for the circumstellar effects that 
create the observed double eclipse pattern, but the stellar photosphere also contributes 
to the brightness variations. This leads to discrepancies (Figure [T]) between the observed 
and simulated photometric light curve, which we are currently investigating. 



4. Spectropolarimetric observations 

We obtained a total of 18 high resolution (R=65000) broadband (370-1040 nm) circular 
polarization spectra of a Ori E. Sixteen of these spectra were obtained in November 2007 
from the Narval spectropolarimeter attached to the 2.2m Bernard Lyot telescope at the 
Pic du Midi Observatory in France. The remaining 2 spectra were obtained in February 
2009 from the spectropolarimeter ESPaDOnS attached to the 3.6-m Canada-France- 
Hawaii Telescope, as part of the Magnetism in Massive Stars (MiMeS) Large Program 
(Wade et al., these proceedings). 

We used the method of Least-Squares Deconvolution (LSD) to derive values of longi- 
tudinal field for each spectrum. LSD describes the stellar spectrum as the convolution 
of a mean Stokes I or V profile, representative of the average shape of the line profile, 
and a line mask, describing the position, strength and magnetic sensitivity of all lines 
in the spectrum. A mask was used for a star with T e ff=23000 K, log g — 4.0, and solar 
abundances except for enhanced helium. From the LSD mean Stokes I and V profiles, we 
calculate the longitudinal magnetic field, B^: 

fr = -2.14x10" (4.1) 
Xgc J[l - I(v)]dv v ; 

(Wade et al. 2000), where A is the weighted average wavelength and g is the weighted 
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average Lande factor in the mask. I c is the continuum value of the intensity profile. The 
integral is evaluated over the full velocity range of the mean profile. 
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Figure 2. New longitudinal field measurements for each individual spectrum of a Ori E (as- 
terisks) along with la error bars. In addition, we have plotted measurements from Bohlender 
et al. (1987), as triangles (hydrogen) and diamonds (helium), and Landstreet & Borra (1978) 
(squares) each with la error bars. 

When compared with historical data from Landstreet & Borra (1978) and Bohlen- 
der et al. (1987), our new measurements agree surprisingly well, indicating stability of 
the global magnetic structure over three decades. The new measurements are also much 
more precise, decreasing the error bars by a factor of four versus the historical data. 
The spectropolarimetric data presented in Figure [2] allow a long baseline to study the 
periodicity of the field as compared with the period derived from photometric studies. 
Using the Scargle periodogram, the entire set of magnetic data gives a period of 1.190842 
± 0.000004 days. The baseline of the data is long enough that the period can be deter- 
mined down to a third of a second. This is strictly an average period over the time frame 
of the observations. As Townsend et al. (2010) have shown, a Ori E is spinning down 
due to magnetic braking at a rate of 77 milli-seconds per year. The associated spindown 
time calculated from observations, 1.34 Myr, agrees remarkably well with the theoretical 
prediction of 1.4 Myr by ud-Doula et al. (2009). 



5. Magnetic Doppler Imaging 

From rotationally modulated line profiles, the surface abundance and magnetic field 
vectors of a star can be mapped using a technique called Magnetic Doppler Imaging 
(MDI), originally developed by Piskunov & Kochukhov (2002) for the Ap star a 2 CVn. 
MDI was adapted from Doppler Imaging (DI), which reconstructs features on the sur- 
faces of stars by inverting a time series of high-resolution spectral line profiles. The code 
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compares observational Stokes parameters with synthetic parameters, using least squares 
minimization. MDI solves the radiative transfer equation to obtain the synthetic line pro- 
files, assuming an unpolarized continuum and LTE. The code can be used for arbitrarily 
complex fields and abundance distributions. 

This paper applies MDI to the Bp star a Ori E. The maps created from circularly 
polarized line profiles show a dipole magnetic field with a polar strength of 9.6 kG, al- 
though higher-pole magnetic field geometries and small scale fields cannot be identified 
without Stokes Q and U information. The surface of the star (as shown in Figure [3j) 
shows spots in carbon and silicon. The helium lines of a Ori E show NLTE effects, espe- 
cially at longer wavelengths. Currently, MDI is fitting helium lines with a partial NLTE 
implementation to include departure coefficients calculated from the NLTE atmosphere 
code TLUSTY (Lanz & Hubeny 2007). It may be possible to fit other elements from a 
Ori E's spectrum, however rapid rotation and lack of available lines affect MDI's ability 
to calculate accurate maps. 




Figure 3. The chemical abundance distribution of a Ori E derived from Stokes I and V profiles of 
the C and SI lines. The star is shown at 5 equidistant rotational phases viewed at the inclination 
angle, i = 75° and vsini = 150 km s _1 . The scale gives abundance as e(Elem) corresponding to 
log(iV Eiem I 'Ntot) • The rotation axis is vertical. 



6. Synthetic photometry 

Krticka et al. (2007) used DI maps of HD 37776 to study the effects of spectroscopic 
spots on the surface of a star on the redistribution of flux. The maps of silicon and 
helium used show that both elements contribute to changes in the continuum of the star's 
emergent flux. In their paper, LTE atmospheres are used to calculate radiative fluxes and 
simulate light curves for each element. The authors note that NLTE effects in lines may 
change the shape of the predicted light curve. Matching the observed photometric light 
curve requires a combination of each elemental light curve. This method has been shown 
to be quite successful at reproducing the light curve of HD 37776, as the star does not 
show an eclipse pattern expected from a magnetosphere. 

For other stars like a Ori E with magnetospheres, Magnetic Doppler Imaging (MDI) 
can be used in the same way to create synthetic light curves for specific elements, taking 
into account the effects of chemical spots on the total flux emitted due to bound-free 
transitions. In §3, we mentioned that the RRM model can provide a synthetic light 
curve for the circumstellar effects on brightness. To fully account for the contributions 
to the emergent flux from a Ori E, we should merge both the light curves synthesized 
from the RRM model and from MDI. Ideally, when we combine both the photospheric 
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and circumstellar brightness variations, we should be able reproduce the details of the 
observed photometric light curve. 

7. Summary 

Massive star magnetospheres are structured, co-rotate with the star, and produce 
rotationally modulated observables. The RRM model of the circumstellar environment 
of a Ori E can reproduce observations relatively well, except for unexplained differences 
in the photometric light curve. From new spectropolarimetric observations we show that 
a Ori E has a strong (9.6 kG), dipole magnetic field that new observations show has 
remained steady in structure over 3 decades. Its longitudinal magnetic field varies from 
-2.3 kG to 2.5 kG. The period derived from these longitudinal field strength variations is 
1.190842 ± 0.000004 days. This precise period is the average over the time baseline, as 
a Ori E's rotation period gains 77 milli-seconds each year due to magnetic braking. We 
can use MDI to map the surface abundance of <j Ori E and synthesize a photospheric 
brightness light curve. By merging the light curves from both the MDI and RRM models, 
we expect to properly account for both the circumstellar and photospheric effects on the 
star's brightness. Using this combined approach, we hope to reproduce the details of the 
observed photometric light curve of a Ori E. 
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